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ABSTRACT 

We use the HB morphology of 48 Galactic GCs to study the radial distributions of the different 
stellar populations known to exist in globular clusters. Assuming that the (extremely) blue 
HB stars correspond to stars enriched in Helium and light elements, we compare the radial 
distributions of stars selected according to colour on the HB to trace the distribution of the 
secondary stellar populations in globular clusters. Unlike other cases, our data show that the 
populations are well mixed in 80% of the cases studied. This provides some constraints on the 
mechanisms proposed to pollute the interstellar medium in young globular clusters. 
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1 INTRODUCTION 


Globular clusters are now known to host multiple stellar popula¬ 
tions, differing in their light element abundances, although iron 
and heavier elements tend to be largely homogeneous (e.g., |Grat-| 
|ton et al^|2012| for a review). It is likely that production of light 
elements is accompanied by enrichment in Helium as well, where 
Na-rich stars from the second generations are more He-rich than 
the Na-poor and O-rich first generation (FG) stars (e.g. |Carretta| 
|et al.|20d7| ). The secondary stellar generations (SG) must have been 
formed during a later starburst, from material polluted by the ejecta 
of massive stars in the original stellar population. Candidate pol¬ 
luters include: intermediate mass AGB stars ( [Izzard et al.| | 2006l 
D^Antona&^entura|2007 1 , fast rotating massive stars ( |Decressm| 
etaL]2007| Krause et al.|2012||2013| >, ejecta from massive binaries 


( |de Mink et al.|2009| > and e^ly disk accretion jBastian et al.|201^ 
where the ejecta from massive stars are accreted onto forming low- 
mass protostars at early times. 

It is generally difficult to distinguish between these scenarios 
purely on the basis of the observed abundance patterns, as they all 
share common characteristics, such as pollution by the products of 
hot bottom burning in massive systems. However, it may be possi¬ 
ble to constrain these models via their effects on the radial distribu¬ 
tion of polluted stars, i.e, the first and later generations. In general. 
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we expect that more enriched stars will tend to reside closer to the 
cluster centres, as the gas needed to fuel star formation tends to sink 
to the cluster core (e.g., |D’Ercole et al.|2008[ l. In the early disk ac¬ 
cretion model, the protostellar disks accrete material as they travel 
through the cluster and therefore one expects a somewhat broader 
distribution. On the other hand, the He-enhanced stars may lose 
more mass during the RGB phase and diffuse outwards because of 
two-body relaxation, yielding a more extended distribution (|Car-| 
|retta et al.|2009| l. The radial distribution of stars in clusters may be 
altered by dynamical evolution. However, simulations by |Vesperini| 
|et al.|p0T^ predict that in many Galactic clusters the second gen¬ 
erations should still be more centrally concentrated than the first 
generation. Opposite, [Decressin et al.|j2008|[20T0) argue that any 
original radial difference is erased after a Hubble time because the 
relaxation times are much shorter than the ages of globular clusters. 

Observationally, [Lardo et alT|j2011^ have used the u — g colour 
to separate red giants belonging to each population (cf. |Milone| 
| et al.| |2008|) and show that the UV-faint, Na-rich second genera¬ 
tion stars are more centrally concentrated. |Carretta et al.| j2009| l 
divide their spectroscopic sample into primordial, intermediate and 
extreme stars based on the degree of enrichment in Na and O abun¬ 
dances and show that the intermediate stars are more centrally con¬ 
centrated, followed by the primordial and then the extreme subpop¬ 
ulations, although |Lardo et^j2011) criticise the ad hoc selection 
of targets. In our previous studies we have shown how selecting red 
giants in the crowded central regions of globular clusters can be dif- 
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ficult and how their photometry may be doubtful, especially from 
small telescopes and in mediocre seeing conditions (|Renzini|1998| 
[Vanderbeke et al.]|2014a|b^. It is perha ps suggestive that the two 
least crowded systems in |Lardo et al.| ( |2011| l are those that show 
no evidence of radial gradients in the distribution of stellar pop¬ 
ulations. Other studies argue that second generation stars may re¬ 
side closer to cluster centres than first generation stars (e^gJSollirm 


et al.|2007]|Carretta et al.|2010a[|Kravtsov et al.|201l] [Johnson & 


Pilachowskij20^, still reflecting the initial segregation. However, 


Dalessandro et al.|([2014[l combined optical and near-UV photom¬ 


etry to study the sub-giant and red giant branches of NGC 6362. 
They concluded that the FG and SG stars share the same radial 
distribution, making it the first system where stars from different 
populations are found to be completely spatially mixed. 

Previous work has used RGB stars as tracers, either from their 
colours or direct spectroscopy (e.g. [Sbordone et al.[[2011^ . Here 
we propose a different approach, based on the distribution of stars 
on the horizontal branch (HB). The dominant parameter determin¬ 
ing the HB morphology is the metal abundance (e.g., [A rp et ak] 
[1952[[Sandage|195^ , while other parameters like age. Helium en¬ 
richment and mass loss, to name a few, may produce anomalously 
red or blue HBs (the ’’second parameter problem” - e.g., [Dotter[ 
[et al.|201^ . Helium abundance variations are linked to the observed 
light element enhancements for the multiple generations in globu¬ 
lar clusters (e.g. Carretta et al.|2006 JDVAntona^^^nturaj2007 


Salads et al.|2008[[Villanova et al.|2009[ [Marino et al.|20H|[Vil 


lanova et al.|2012||Monelli et al.|20T3l [Salads et al.|2013|[Mi 


et al.[[20f4l [Gratton et al.[[2014l [Mucciarelli et al.[[2014| [Gratton 


et al.[[20f5) [Milone[[2015[ > and may be taken to trace abundance 
variations in CNONa, even if these by themselves may not affect 
the morphology of the HB. 

[Carretta et al.H2010b[ l show that the maximum temperature of 
the HB correlates with the range of Na/0 abundance in red giants. 
[D’Antona et al.| ( |2005^ and [lannicola et ar] ( [2009[ l compare their 
data on NGC 2808 with synthetic models and demonstrate how the 
He abundance increases bluewards along the HB in their colour- 
magnitude diagrams (CMDs). The models of [ Joo & Lee[ ( [20 1 3} also 
link the HB morphology directly to He enrichment and the presence 
of multiple stellar generations, although this is a 2nd order effect on 
the dominant metallicity which is the 1st parameter affecting the 
HB. [Marino et al.[ ( [2011[ l find that in M4 the red HB (RHB) stars 
are O-rich and Na-poor, while the blue HB (BHB) stars are O-poor 
and Na-rich; the abundance patterns suggest that the elements were 
produced during hot bottom burning via the CNO cycle and its high 
temperature NeNa and MgAl branchings (e.g., |Clayton|1968^ and 
the resulting stars will also be Helium enhanced. Therefore, the 
distribution of stars on the blue HB within each cluster may be used 
as a proxy for the distribution of stars belonging to each generation. 

We exploit the relation between light element enrichment and 
HB morphology, mediated by He abundance, to study the spatial 
distribution of stellar populations selected on the HB. We essen¬ 
tially take temperature ’cuts’ along the HB and argue that these cor¬ 
respond to increasing contributions from He-enriched secondary 
generations of stars. For example, [lannicola et al.[j2009) can sepa¬ 
rate the HB of NGC 2808 into three groups, each assumed to cor¬ 
respond to the three (or more) main sequences observed by [Pi-[ 
[otto et ar] ( [2007^ . Observations in near ultraviolet generally show 
numerous gaps on HBs ( [lannicola et al.[[2009^ , even if these are 
continuous in the optical, and these may be used to separate the 
multiple stellar generations. [lannicola et al^ ( [2009^ selected about 
2000 HB stars based on both ground-based and HST observations 
of NGC 2808. Relative fractions of cool, hot and extreme HB stars 


do not change radically along the radial profile of the cluster. There¬ 
fore, this is evidence against the presence of radial differences be¬ 
tween stellar subpopulations with different He abundances. [Kund^ 
[et al.[ ( |20T3t studied the cumulative fraction of extreme HB stars in 
M 22 (NGC 6656) but could not draw strong conclusions regarding 
differences in the radial distributions of blue and extremely blue 
HB (EHB) stars. [Gratton et aL] ( [2014^ studied M 22’s HB spectro¬ 
scopically and find some suggestive evidence that SG stars are more 
concentrated than FG stars. However, their study does not cover the 
extreme HB stars, due to the restrictions on the temperature range. 

Here we use our Galactic Globular Cluster Catalog (G2C2) 
photometry ( [Vanderbeke et al.|2014a|b[ see Vanderbeke et al. 2015, 
in preparation, for a discussion of the colour-magnitude diagrams 
of all clusters in the G2C2 sample) to address this issue. Our large 
and homogeneous photometry is sufficiently deep that we can trace 
the HB to the level of the main sequence turnoff in almost all clus¬ 
ters and allow us to study the radial distributions of stars at dis¬ 
tances beyond several core radii from the cluster centre in many 
instances (the CTIO FOV is 13.6'on the side and has a resolution 
of 0.396”per pixel). The main goal of this study is to use the distri¬ 
bution of stars on the HB as a proxy for the distribution of stellar 
populations and use these stars to trace the radial density profile 
of the various stellar generations. We wish to test whether second 
generation stars are more centrally concentrated than the more pri¬ 
mordial objects and if possible consider the effects of dynamical 
evolution over the past Hubble time. Our homogeneous and wide- 
field photometry proves ideal for this task. 


2 SAMPLE SELECTION AND METHODOLOGY 


The basic data for this paper derive from our homogeneous pho¬ 
tometry of Galactic globular clusters in Papers I and II (|Vander-| 
[beke et al.|20I4a|bl (. The full colour-magnitude diagrams for these 
clusters will be fully discussed in a subsequent publication. 

Because we are interested in identifying objects with extended 
HBs, and using colour cuts on the HB to select stars in temperature 
ranges so that each group contains stars from the primordial or sec¬ 
ondary generations, we focus on clusters whose HBs have blueward 
extensions (see e.g. [Mackey & van den Bergh|2005[[Lee et al.|2007[ 
[Potter et al.|2010) . In Fig. [T] we present histograms of the HB in¬ 
dices of the GCs studied in this paper. We included all GCs where 
the HB index is larger than 0.7, whenever good quality data (deep 
CMDs, low contamination) was available. 

While blue stragglers may contaminate counts of extreme HB 
stars, in most cases only the hottest and more massive blue strag¬ 
glers would be confused with HB stars, and these are generally very 
few in number. The location of the field stars in the CMD varies 
from cluster to cluster. However, the coinciding HB region would 
appear less centrally concentrated by the field star contamination. 

Some clusters with rather red HB morphologies were also in¬ 
cluded, because these GCs either belong to a second parameter pair 
(e.g. NGC 362, NGC 6171), or there is existing evidence of mul¬ 
tiple populations with differing radial distributions (e.g. NGC 104, 
NGC 6362), or visual inspection revealed some BHB/EHB stars 
though the HB is strongly dominated by the RHB (e.g. NGC 1261). 
After correcting our CMDs for foreground extinction using 
Cardelli et al. ^1989^ reddening law (using E{B — V) val- 


the 


ues from the 2010 version of [Harris|1996t , we subdivide the HB 
in ranges, according to colour, where the first or second generation 
stars are expected to prevail (cf., [lannicola et al.|2009|[Marino et ^ 
[201 1[|. This is unfortunately not straightforward when using optical 
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colours. [Joo & Le^j2013^ predict different colours for the differ¬ 
ent generations residing in M 22 (NGC6656) and NGC1851 and 
show that the generations can have partly overlapping colour and 
magnitude ranges. Therefore, every possible cut will only be an 
approximate separation of the populations. In this study we make 
the assumption that the majority of RHB/BHB stars represent the 
stars with primordial abundances, while the majority of EHB stars 
correspond to the enriched population. Either sample is of course 
contaminated to some extent. 

For the sake of homogeneity, we define colour cuts which we 
will use systematically for all clusters using NGC 1851 to illustrate 
the different steps in our analysis (see Fig.|^. This cluster is known 
to host multiple stellar generations and is famous for its unusual 
bimodal HB morphology, double RGB and SGB, possible [Fe/H] 
spread, CN bimodality and variations in light and s-process element 
abundances ( [Saviane et al.|19^ | Walker|1998| |Jo o & Lee|2013| >. 
ICarretta et al. ( |201 1| ) demonstrate that the metal-poor stars in their 
sample are more concentrated than the metal-rich stars. 

Now we define colour cuts for the different HB subgroups. 
|Dalessandro et all]j2011^ showed the temperature distribution of 
NGC 2808 HB stars and found an effective temperature of about 
7000 K as the limit between RHB and BHB stars, close to the blue 
edge of the RR Lyrae instability strip jMoni Bidin et al.|2012^ . This 
temperature corresponds to a black body colour of g—z ~ 0.29 and 
was used as a vertical cut to separate RHB and BHB stars. The gap 
between BHB and EHB stars is located at about 20000K ( |Dalessam| 
|dro et al.|201 l[[Moni Bidin et al.|2012y However, this Te// corre¬ 
sponds to a black body colour of g — z ^ —0.9, which is bluer than 
any of the HB stars in our CMDs. [Dalessandro et al.| ( [2011[ > demon¬ 
strated that the EHB temperatures derived from a combination of 
optical filters only can be underestimated by lOOOOK and more (be¬ 
cause HB sequences with different initial He mass fraction overlap 
for Te// > lOOOOiT). Therefore, we are forced to introduce an ar¬ 
bitrary colour cut below the ’bend’ of the HB at g — z ~ —0.4, cor¬ 
responding to an effective temperature of 1 lOOOK. Because the HB 
becomes vertical at blue colours, we do not make a pure ’’vertical” 
colour cut for the EHB stars. We fit a second degree polynomial to 
the HB stars and determine the line perpendicular to this polyno- 



Figure 3. Zooming in on the HB of NGC1851. The solid line is a quadratic 
fit to the HB stars. The dashed lines are the colour cuts applied to separate 
RHB (cooler/redder than g — z = 0.29, corresponding to ~ 7000iT), BHB 
and EHB (hotter/bluer than the dashed line passing through g — z = —0.4 
(corresponding to ~ llOOOiT), which is perpendicular to the solid line) 
stars. 


mial passing through g — 2 ~ —0.4. In Fig. we zoom in on the 
HB region of NGC 1851. We illustrate the fitting procedure by the 
solid line and the different cuts by the dashed lines. 

Our CMDs are derived from data observed with small ground- 
based telescopes (both Sloan Digital Sky Survey (SDSS) and Cerro 
Tololo Inter-American Observatory (CTIO) 0.9m telescope), as 
presented in |Vanderbeke et ar|j2014a| l. Therefore, these suffer from 
incompleteness, especially close to the cluster centre. We show the 
number of stars as a function of distance from the cluster centre in 
0.5' annuli as a function of luminosity in Fig.|^ The HB spans a 
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Figure 2. NGC1851: Left panel: a gz CMD of NGC 1851 from our data, corrected for foreground extinction (Cardelli et al.|1989| using E{B — V) from 
the 2010 version of[H arris The errors shown are the median photometric errors. These errors are the errors given by DAOPHOT jStetson|1987) and 

do not include the photometric calibration errors or the systematic photometric uncertainty found in Paper I (as these would dominate the error budget). We 
systematically exclude stars with a magnitude error larger than 1. Right panel: cumulative fractions of the radial distributions of the HB regions. The number 
of stars within the inner and outer limiting completeness radii are given for the different HB regions. 



Annuli [arcmln] 


+ Q I I 

O 1 4- 1 5 mag 
1 5— 1 6 mag 
A; 1 6 - 1 7 mag 
O 17-18 mag 
X 18-19 mag 

♦ HB 15-16 mag 
HB 16-17 mag 

♦ HB 17-18 mag 

♦ HB 18-19 mag 


Figure 4. NGC1851: radial distribution of the number of stars contained in 0.5'annuli. The half-light radius is indicated with dashed line, the dotted line 
represents the radius where pai1 of the annulus is outside the CTIO field of view. 
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magnitude range 15.5 < 2 < 18, so we need similar completeness 
levels for that magnitude range. Therefore, we choose to consider 
only the HB stars beyond an inner limiting radius of 1.5'. In Fig.|^ 
the half-light radius is indicated with dashed line, the dotted 
line represents the radius where part of the annulus is outside the 
13.6' CTIO field of view. Note that the inner radial limit is well be¬ 
yond the half-light radius. As the outer limiting radius, we choose 
8' (hence including almost the entire CTIO field of view). 

In general, we use all stars to determine the completeness at 
each radius. However, the total number of stars is dominated by the 
numerous RGB stars, while the focus of our study is on the HB 
stars, which are much bluer. Therefore, we still need to check that 
the total number of stars is representative to make the completeness 
cuts for the HB, so, that no bias for blue stars was introduced by the 
characteristics of the CCD. The filled symbols in Fig. [^represent 
the HB stars within the given magnitude ranges. HB stars of differ¬ 
ent magnitude ranges within the previously determined inner and 
outer limits show similar completeness levels (although they suffer 
from low number statistics). Similar conclusions were drawn for 
another dozen clusters spanning a variety of globular cluster prop¬ 
erties (including mass, HB morphology, etc.). Therefore, we con¬ 
clude that we can safely use all the stars, even if dominated by the 
RGB stars, to determine the completeness as a function of radius. 

We define and colour-code the different HB regions in the left 
panel of Fig. and compare the cumulative radial distributions of 
the stars between the inner and outer completeness radii in the right 
panel of the same figure. Similar figures for all other GCs can be 
found in the online appendix. 

We perform a two-sided KS test to compare the radial distri¬ 
butions of the HB stars as defined above (RHB, BHB and EHB). 
The test returns no significant difference between the radial distri¬ 
butions of RHB, BHB and EHB stars. In Table[T]we show an extract 
of the KS statistics to guide the reader. The complete table can be 
found in the online appendix. Note that our sample also includes 
u) Cen (NGC 5139), an object that it is much more complex than a 
typical mono-metallic globular cluster (e.g., |Villanova et al.|2014^ . 
A more detailed analysis of effects of the [Fe/H] spread on the HB 
selection criteria is beyond the scope of this study. 

[Saviane et al.| ( |199^ also studied NGC 1851 and found evi¬ 
dence that the radial distribution of the blue HB stars is different 
from that of the red HB and sub-giant branch stars. Their Fig. 11 
shows that blue horizontal branch stars are more centrally concen¬ 
trated than red horizontal branch stars, while our figure suggests 
that only the EHB stars are slightly more concentrated than the 
RHB stars. Our KS statistics also indicate that the difference in 
cumulative radial distributions is not very significant.|Saviane et al.| 
( |1998^ also presented HST imagery for the inner 25 arc seconds. 
These data did not suggest any significant radial difference for the 
blue and red HB stars. [Milone et al.| ( |2009t also studied NGC 1851 
and did not find any radial stellar population gradients, in agree¬ 
ment with our results. 


3 RESULTS 

In Eig. we show the distribution of the KS probabilities (from 
CTIO and SDSS data). The distribution of the KS probabilities rel¬ 
ative to the EHB is given by the dashed histogram. The distribution 
is trimodal: in some clusters there are strong differences in the ra¬ 
dial distribution of HB stars within each colour range; others show 
either mild or no differences. 

More than 80% of the clusters do not show evidence for dif- 



Figure 5. Histogram of the KS probabilities. The histogram without dashes 
presents the distribution of all clusters, while the dashed histogram repre¬ 
sents the subsample for which we consider the EHB. See Tablefor the 
specific values. 

ferent radial distributions along their horizontal branch. This sug¬ 
gests that the different stellar populations have similar radial dis¬ 
tributions once on the HB, at least for the considered radii, in con¬ 
trast to the study by |Lardo et al.| ( |2011^ who probed the red gi¬ 
ant branch stars. Moreover, we find no significant radial distribu¬ 
tion difference for several clusters for which |Lardo et al.| j2011| > 
found differences in the radial distributions of first and second gen¬ 
eration stars: NGC 5024, NGC 5904, NGC 6205, NGC 6341 and 
NGC 7089 (M 2). If the position of stars on the HB is related to 
their He abundance (and indirectly to the CNONa anomalies), this 
result is puzzling. In agreement with |Lardo et al.|f2011^ we find no 
difference in the radial distributions of stars for NGC 5466, while 
in NGC 5272 we only find a small difference in radial distribu¬ 
tions (Probics ~ 0.2 for the EHB comparisons), while |Lardo et al.| 
( |2011^ found a significant difference. For NGC 7078 we find a dif¬ 
ference in the radial distributions for both the RHB-EHB and RHB- 
BHB comparisons as do |Lardo et al.|p011| >. For NGC2808 we do 
not find a radial difference between RHB and EHB, in agreement 
with [lannicola et aL]j2009^ . Eor NGC 6362, we found that RHB 
and BHB stars have the same radial distribution, agreeing with |Da-| 
[lessandro et al.| ( (2014| l. However, the EHB stars are somehow more 
concentrated than the RGB stars (with Probirs ~ 0.13). 

NGC 288, NGC 362 and NGC 6218 are particularly inter¬ 
esting objects (all with [Fe/H\ ~ —1.3), as these are the first 
known systems for which the second generation appears less con¬ 
centrated than the first generation, based on our data. Photometric 
or spectroscopic follow-up studies are needed to confirm or dis¬ 
prove these radial distributions. Disc shocking may be a main con¬ 
tributor to the peculiar radial distributions of the different popula¬ 
tions in NGC 288 ( [Kruijssen & Mieske|2009^ . 

3.1 A pure red HB cluster: NGC 104 

IGratton et al.| ^2013^ demonstrated that only the reddest HB stars 
in NGC 104 (47 Tuc) can be considered as FG stars, the bluest 
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Table 1. Extract of the KS statistics table for radial distributions of different HB regions. HBi and HB 2 denote the considered HB regions, as defined in the 
linked figure. Ni and N 2 give the number of the stars in both HB regions. g gives the Kolmogorov-Smimov statistic and Probes presents the significance 
level of the KS statistic. Small values show that the cumulative radial distribution of HBi stars is significantly different from HB 2 stars. The last column 
indicates if the CMD is based on CTIO or SDSS data. 


ID 

HBi 

HB 2 

Ni 

N 2 

Dks 

Probxs 

Fig. 

NGC 104 

red 

blue 

245 

230 

0.198 

0.000 


1 

CTIO 

NGC288 

BHB 

EHB 

17 

67 

0.345 

0.061 



CTIO 

NGC362 

BHB 

EHB 

10 

26 

0.223 

0.814 



CTIO 

NGC362 

RHB 

EHB 

67 

26 

0.372 

0.008 



CTIO 

NGC362 

RHB 

BHB 

67 

10 

0.301 

0.344 



CTIO 

NGC 1261 

BHB 

EHB 

9 

6 

0.556 

0.140 



CTIO 

NGC 1261 

RHB 

EHB 

93 

6 

0.226 

0.897 



CTIO 

NGC 1261 

RHB 

BHB 

93 

9 

0.459 

0.043 



CTIO 

NGC 1851 

BHB 

EHB 

27 

18 

0.222 

0.603 


2 

CTIO 

NGC 1851 

RHB 

EHB 

57 

18 

0.202 

0.583 


3 

CTIO 

NGC 1851 

RHB 

BHB 

57 

27 

0.183 

0.526 


3 

CTIO 


ones are enriched in Na and depleted in O. In Fig. we show the 
CMD and radial cumulative fractions for NGC 104. Although our 
’’standard” approach is tailored to clusters with extended HBs, it is 
possible to make an arbitrary cut in 47 Tuc and study the cumula¬ 
tive radial distributions. We recover the |Gratton et al.| ( |2013t result 
and find that the blue part of the RHB is statistically more centrally 
concentrated than the red part (with probKS = 0.0001). It further 
confirms the results of |Nataf et al.H20fT} , who found evidence for 
a centrally concentrated, He-rich SG. However it may be an ad hoc 
result and one should apply this method to metal-rich or red HB 
clusters with caution, although the success of this experiment sug¬ 
gests that our approach does select ranges where one or the other 
population is more significant. 


4 SUMMARY 

Our data show that the distribution of stellar populations is more 
complex than expected. While most GC formation scenarios pre¬ 
dict that the second and further generations will have different ra¬ 
dial distribution from the primordial cluster populations, usually in 
the sense of being more centrally concentrated, we find that there 
is little evidence that this is generally true, if we use HB stars as 
tracers of He enriched population, a feature that should accompany 
the light element enhancement typical of the secondary stellar gen¬ 
erations. 

Although in some cases we expect that the stellar populations 
will have been thoroughly mixed by dynamical evolution, including 
the possible effects of disk shocking, simulations of |Vesperini et al.| 
( |2013^ predict that original population gradients will not have been 
erased by the present time in many GCs. 

In that respect, NGC 288, NGC 362 and NGC 6218 are of 
specific interest. For these systems, the HB stars linked with the 
second generation are significantly less concentrated than the RHB 
stars. Photometric or spectroscopic follow-up studies are needed to 
confirm or disprove these radial distributions. 

One caveat on the present work is that we have to assume that 
the spread in colour on the HB in each cluster is due to variations 
in helium abundance and can be related to the various stellar gener¬ 
ations known to be present in these objects. While this seems to be 
a reasonable assumption there are clusters where no helium vari¬ 


ation is apparent on the HB, even though chemical anomalies and 
multiple populations are present. 

Taken at face value, our results are not fully consistent with 
current enrichment scenarios dominated by AGB stars or fast ro¬ 
tating massive stars, as these would generally produce more highly 
concentrated second generation stars. Nevertheless, this points to 
the necessity of improving our theoretical understanding and mod¬ 
elling of multiple stellar populations in clusters, as well as defining 
a consensus tracer population. Hubble Space Telescope photome¬ 
try (e.g. |Piotto et al.|2015| ) allows to extend our analysis into the 
cluster centres and could potentially provide additional clues. 
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